Spin polarization of half-quantum vortex in systems with equal spin pairing. 
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We present a variational analysis for a half-quantum vortex (HQV) in the equal-spin-pairing 
superfluid state which, under suitable conditions, is believed to be realized in Sr2Ru04 and '^He-A. 
Our approach is based on a description of the HQV in terms of a BCS-Hke wave function with a 
spin-dependent boost. We predict a novel feature: the HQV, if stable, should be accompanied by 
a non-zero spin polarization. Such a spin polarization would exist in addition to the one induced 
by the Zeeman coupling to the external field and hence may serve as an indicator in experimental 
search for HQV. 



When '^He is liquified and cooled into the millidegree 
regime it enters a new phase which has been proved to 
possess a spin triplet paired condensate. There also ex- 
ists by now a growing body of experimental evidence that 
Sr2Ru04 below 1.5 K is a spin triplet superconductor [l[. 
This implies the possibility of many interesting phenom- 
ena not expected in systems with spin singlet pairing. 
One of them is the existence of half-quantum vortices 
(HQV's) in the equal-spin-pairing (ESP) state of the spin 
triplet condensate 0, 0] . 

The pairing symmetry for ^He is well established and 
the so-called A phase is confidently believed to realize an 
ESP spin triplet state [1, While there is, to the best 
of our knowledge, no unambiguous observation of HQV 
in '^He-A there is a strong theoretical argument in 
favor of their existence at least under some assumptions 
on the geometry of the experiment. On the contrary, the 
pairing state of Sr2Ru04 is currently poorly understood 
and the observation (or not) of HQV, along with other 
experimental information, would facilitate identification 
of the underlying pairing symmetry. In addition, there 
is a significant interest in HQV for topological quantum 
computing [l3| . 

It should be noted right away that the identification 
of the superconducting phase of Sr2Ru04 with an ESP 
state does not by itself guarantee thermodynamic stabil- 
ity of HQV in this compound. Even under the assump- 
tion of negligible spin-orbit coupling, the kinetic energy 
of unscreened spin currents which accompany HQV dis- 
favors its formation vis-a-vis the formation of a regular 
vortex where electromagnetic currents arc screened over 
the length of the penetration depth. Such an unfavorable 
energy balance can be avoided by limiting the sample size 
to a few microns 0. This however further complicates 
the experimental detection of HQV in Sr2Ru04. 

One of the most direct ways for detection of HQV is 
to look for spin currents which circulate around it. The 
usual techniques for spin current detection arc based on 
the accumulation of spin and their straightforward ap- 
plication to this situation seems to be difficult. One can 
however use a fact that spin currents generate electric 
field. A very rough conservative estimate shows that for 
a ring of size 1/i the quadrupole electric field generated 



by the spin currents of HQV will create a potential differ- 
ence of InV across the ring ~ quite small, but not beyond 
the capabilities of current experimental techniques. 

There is also a possibility of detecting HQV by looking 
for specific features in the magnetization curves of small 
rings made of Sr2Ru04. Experiments in this direction 
are currently underway in the Budakian group at UIUC. 

In this paper we suggest an apparently new effect 
which may be utilized for the detection of HQV both in 
Sr2Ru04 and in ^He-A. The effect consists in the pres- 
ence of an effective Zeeman field in the HQV state of the 
ESP condensate. In thermodynamic equilibrium such an 
effective Zeeman field will produce a non-zero spin po- 
larization in addition to that created by external fields. 
In particular, such a spin polarization would exist even 
in the absence of external Zeeman coupling provided the 
condensate is in the HQV state. At the same time this 
field would not exist in a normal vortex state thus allow- 
ing one to distinguish between the two. For a ly^t ring 
of Sr2Ru04 the magnitude of the eflFective Zeeman field 
is about 10 Gauss and, taking the spin susceptibility to 
be of order 10~^emu/mol [l[, the spin polarization pro- 
duced by such a field can be seen in Ti or even Knight 
shift measurements. 

We start by noticing that in the ESP state the Cooper 
pair is always in a linear superposition of states in which 
both spins in the pair are either aligned ("up") or an- 
tialigned ("down") with a given direction in space. The 
corresponding many-body wave function for a system 
with N /2 pairs which are condensed into the same two- 
particle state characterized by functions ip'\ and Lpi can 
be written as 

*ESP [<^T(i-i,r2)|Tt) + '^i(i-i,r2)|ii)] ... 

|_(pi(r7v-i,r7v) ITT) +</3x(rN-i,rjv) |ii.)J j, 

where A is the antisymmetrization operator with respect 
to particles' coordinates and spins. In the weak cou- 
pling limit, provided the pairing interaction conserves 
spin, the up and the down spin particles can be consid- 
ered as independent subsystems. In this case the HQV 
state of the ESP condensate has a simple physical inter- 
pretation: It is a state in which the two spin systems 
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have different winding numbers, i.e. accommodate differ- 
ent number of vortices. 

To avoid comphcations related to the presence of the 
vortex core we specialize to an annular geometry. Let R 
be the radius of the annulus and d be the wall thickness; 
it will be assumed that d/i? <C 1 so that effects of order 
d/R or higher can be ignored. Then specializing to the 
zero-temperature case and choosing the spin axis along 
the symmetry axis of the annulus, a conceptually simple 
ansatz for the HQV state of the condensate is 

*HQV - exp { ^ 5] 0, + ^ 5] e,} vl-ESP, (2) 

where 6i denotes the azimuthal coordinate of the i-th 
particle on the annulus. The integer is a projection 
of the angular momentum of the a-th component of the 
wave function of the pair on the symmetry axis; the coef- 
ficient 1/2 in the exponent reflects the fact that is the 
momentum of the pair, and, the case £| = corresponds 
to a regular full vortex. As can be seen from the above, 
state ^'hqv is obtained from the initial state ^'esp by 
a uniform spin-dependent boost. While there might be 
doubts that the actual HQV is described by such a sim- 
ple form, it nevertheless should be considered as a good 
starting point for a variational analysis. 

In the d-vector formalism ^hqv s^s written above pro- 
duces a d vector which lies in the plane perpendicular 
to the spin axis i.e. to the symmetry axis of the annu- 
lus. For an annulus made of single crystal Sr2Ru04 with 
the c-axis along the symmetry axis, this corresponds to d 
being in the a6-planc of the crystal. Although this config- 
uration is not favored by the spin-orbit interaction, there 
arc theoretical indications that even a very small exter- 
nal magnetic field along c-axis can stabilize it Q . Similar 
considerations also apply to "^He-A, and in what follows 
the in-plane position of the d vector will be assumed. It 
should be emphasized however that our qualitative con- 
clusion about non-zero spin polarization in HQV does 
not depend on this assumption. 

Knowledge of the state ([2]) allows one to obtain a vari- 
ational energy of HQV which can then be minimized to 
yield its detailed structure. For definitencss from now on 
we will consider only charged systems. In this case the 
appropriate for minimization thermodynamic potential is 
Gibbs energy: 



(3) 



where (Ti.) is the expectation value of the Hamiltonian 
of the system and B and H are the magnetic field and 
induction respectively. For an annulus in a shape of an 
infinite cylinder with the fields along the symmetry axis 
H is the external magnetic field and B is the field inside 
the cylinder. 

As it will be seen below the actual form of Ti is crucial 
for the stability of the HQV. In the simplest case Ti, can 



be taken to contain only the reduced BCS Hamiltonian 
'^BCS with the spin triplet pairing term. However such a 
choice of Ti combined with the ansatz ([2]) never makes 
HQV thermodynamically stable; at best the HQV, in 
which ^1 7^ £i, is degenerate with a full vortex £| = 
£l at the transition point between states with different 
vorticities. To lower the energy of HQV below that of a 
full vortex one needs to account for strong interparticle 
forces. This can be done in the framework of Fermi liquid 
theory which is also applicable in the superconducting 
state With that purpose we write the Hamiltonian 
of the system as 



H — Hbcs + "Hph- 



(4) 



Here T^fl describes energy corrections due to Fermi liq- 
uid effects and Wbcs is a reduced BCS Hamiltonian with 
spin triplet pairing term representing the weak coupling 
part of the theory. We will first evaluate the expectation 
value of the weak coupling Hamiltonian on the state ([2]). 
It can be written as a sum of three terms which have 
different physical origins: 



^-Bcs — Ef) + Es + T. 



(5) 



The first term in the equation above is the energy con- 
tribution coming from the internal degrees of freedom 
of Cooper pairs. For the radius of the annulus R much 
larger than the BCS coherence length this contribu- 
tion will depend on neither the center of mass motion of 
the Cooper pairs, i.e. on quantum numbers £| and i^, nor 
the magnitude of the magnetic field . Assuming that 
we are dealing with a big enough annulus this term will 
not be included in the subsequent considerations. 

The second term is the spin polarization energy of the 
system. Let N„ be the number of particles with spin 
projection a. Defining S* as a projection of the total 
number spin polarization on the symmetry axis 



S=iN^~N^)/2, 



(6) 



and gs as the gyromagnetic ratio for the particles in ques- 
tion, the spin polarization energy takes the following form 



Es 



2xesp 



- gsfiB B ■ S, 



(7) 



where xesp is the spin susceptibility of the ESP state 
calculated in the weak-coupling limit. It should be noted 
that at this point the total spin polarization 5 is a vari- 
ational parameter with the actual value of S to be found 
by the minimization of energy. 

The third term on the r.h.s of cqn. ^ is the kinetic 
energy of the currents circulating in the system. Let $ 
be the total flux through the annulus and $o = hc/2e be 
the flux quantum. Introducing the notation 



$/$o, 



it 



(8) 
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one obtains for T the following expression 



T = 



8m*i?2 



{(4'* + 4'p)^ + 



(9) 



where m* is the effective mass of the particles due to 
Fermi liquid corrections [l^ . In the expression above 
the first term in the brackets is proportional to the to- 
tal number of particles N = + and is thus fixed 
for given values of £s, igp- The second term is propor- 
tional to the spin polarization S and creates an effective 
Zeeman field in the HQV state due to a mismatch be- 
tween velocities of the up and down spin components. 
The value of this field and hence the magnitude of the 
thermal equilibrium spin polarization should be found by 
energy minimization. However, as have already been em- 
phasized, minimization of Ebcs (or corresponding Gibbs 
potential G when self inductance is important) does not 
produce a stable HQV; at best the HQV is degenerate 
with a full vortex at the transition point at which the 
effective Zeeman field vanishes due to vanishing of £s*- 

To make an HQV stable one needs to go beyond the 
weak coupling Hamiltonian and introduce strong cou- 
pling effects. This can be done in the framework of Fermi 
liquid theory, in the way indicated by eqn. For that 
we need to calculate the change of the Fermi liquid en- 
ergy E'fl caused by the presence of spin and momentum 
currents in the HQV state. These currents are generated 
by the spin-dependent boost ([2]) and can be expressed 
in terms of spin-up and spin-down quasiparticle distri- 
butions. Using the standard formalism of Fermi liquid 
theory one obtains: 



44$4p(^^i 



(10) 



Here {dn/de) is the density of states at the Fermi level 
and Zq , Zi and Fi are Landau parameters [3| ■ The first 
term, proportional to Zq, is the energy cost produced by 
a spin polarization and the rest describes Fermi liquid 
corrections due to the presence of the currents. 

It is worth pointing out that expressions Q and pD)) . 
which describe energy transformation under the spin de- 
pendent boost ([2]), can also be written down in terms of 
momentum and spin currents and are limiting forms of 
more general transformation rules given in, e.g. 



Now we are in a position to find the equilibrium spin 
polarization in the HQV state. Minimizing the energy 
Ebcs + -^'FL with respect to 5* we obtain that in cc^uilib- 
rium 



(11) 



where x is the spin susceptibility of the system which, 
up to terms of order e^^h? /2m* , is the spin suscepti- 
bility of the ESP state with Fermi liquid corrections; for 



^He-A at low temperatures the value of x is about 0.37 
of the normal state susceptibility The other quan- 
tity of interest in eqn. pT|) . the Zeeman field B, has two 
contributions: 

B = B + B,s, (12) 

which are the external Zeeman field B and the effective 
Zeeman field i?off caused by the presence of spin currents: 



Br 



${l + Fi/3 + Zi/12}. (13) 



2TO*i?2 

In thermal equilibrium the effective field is a periodic 
function of the total flux $ with period $0 and changes 
its sign at flux values equal to half-integer number of flux 
quanta. Since at least some of the constants entering 
eqn. (|13p are currently not known for Sr2Ru04 it is not 
possible to give an accurate prediction of the field's mag- 
nitude. It is, however, of order ii^^h?/2mR'^ = $o/7ri?^; 
since the first HQV, if stable, exists at about the same 
value of the external field, this means that the spin po- 
larization produced by the effective field is comparable to 
that induced by the external field. It is this phenomenon 
which may provide additional ways for the experimen- 
tal detection of HQV. Its signature would be a sawtooth 
contribution given by eqn. p3p . to the otherwise linear 
field dependence of the Zeeman spin polarization. 

Taking into account both types of spin polarization 
and omitting the internal energy contribution (cf. the 
discussion after eqn. ([5])), the energy of the system E = 
(Ti.) can be written as 



E 



Zi/12 



l + Fi/3 



(14) 



where m is the bare particle mass related to m* by the 
usual relation of Fermi liquid theory. For reasonable val- 
ues of the external field the contribution of the spin po- 
larization energy given by the first term on the r.h.s. 
of eqn. (fT4|) relative to the total energy E is of order 
fv^e^^ /2mR^ and thus can be safely ignored for the anal- 
ysis of the stability of HQV. 

The region of stability of the HQV depends on (1 + 
Zi/12)/(1 -I- Fi/3) which is a zero temperature value for 
the ratio of spin superfluid and superfluid densities Psp/Ps 
[J. The stability criterion found by direct minimization 
of (fH)) yields the condition Psp/Ps < 1 [l3|, which is usu- 
ally fulflUed. However, as has been pointed out by Chung 
et al. the self inductance effect, whose treatment ne- 
cessitates the use of the Gibbs potential ^ constructed 
out of the energy ([M]). replaces this condition with a 
much more stringent one. In particular, for a cylindrical 
annulus the stability of HQV requires that 

Psp/Ps<{l + Rdl2\iy\ (15) 

with Al denoting the London penetration depth. The 
value of Psp/Ps in Sr2Ru04 is currently unknown, how- 
ever condition (fT51) makes the existence of HQV in large 
rings practically impossible. 
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The physical interpretation of the stabihty condition 
(fTS)) is quite transparent: While the electromagnetic 
currents which accompany both the full and the half- 
quantum vortex are well screened in the bulk, the spin 
currents which are present only in the HQV are not, pro- 
ducing an additional energy cost over the full vortex. To 
mitigate such a cost and make HQV stable one needs to 
reduce the spin current energy by reducing either spin 
superfluid density psp or the "effective volume" Rd/2\^ 
over which the spin currents flow such that the condition 
(fT^ is satisfied. 

It is further to be remarked that in an annular geome- 
try where one does not have to deal with the vortex core, 
the stability condition Psp/ Ps < 1 is still valid for a neu- 
tral ESP superfluid such as '^Hc-A. In '^He-A the ratio 
Ps-pl Pa is known to be well below 1 for all temperatures 
below critical, hence making the existence of HQV pos- 
sible in a large part of the phase diagram. By contrast, 
recent numerical results |12l | in a solid cylindrical geom- 
etry claim that HQV exists only in a high field region 
and at temperatures sufficiently close to the transition. 
Wc believe, however, that the narrowness of the region 
of HQV stability obtained in [l^] is due to the omission 
of the Fermi liquid effects from the consideration. 

The real-life question about the thermodynamic sta- 
bility of HQV is complicated and may depend on many 
factors not included in the preceding discussion. Among 
these are deviations from the annular geometry, inclusion 
of spin-orbit interaction and the possibility of the d vec- 
tor lying in the plane other than ah. It should however be 
emphasized that once the HQV has been stabilized we do 
not expect our qualitative conclusion about the presence 
of the spin polarization to be altered by the aforemen- 
tioned factors since this conclusion originates from one 
of the defining properties of the HQV, namely velocity 
mismatch between different spin components. This ve- 
locity mismatch shifts the chemical potentials of up- and 
down- spin components by an amount of order h^/niR^ 
which, in thermal equilibrium, produces an effective spin 
polarization. 

In conclusion, we have shown that the thermal equi- 
librium state of the half-quantum vortex in the annular 
geometry should be spin polarized. This effective spin 
polarization is a periodic function of flux and contributes 
additively to the spin polarization induced by the exter- 



nal Zeeman coupling. The magnitudes of the two con- 
tributions are comparable, thus making the effect po- 
tentially observable through, e.g., NMR measurements. 
This suggests a new way for the experimental detection 
of half-quantum vortices. 
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